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Introduction 

Let / be a homogeneous ideal in A := k[Xo, . . . , X n ] (where k is a field) given in terms 
of its generators 

where fa is a form of degree di. 
This / defines a scheme 

Zj C P n (fc), 

and there is a one to one correspondance between the subschemes of P n (k) and the ho- 
mogeneous ideals up to saturation (the saturation /* of I consists of elements / in A such 
that for some m, X™ f is in / for all i). 

We will mention few ideas for algorithms to compute geometric informations on Zi 
from the generators of / and speak about one aspect of their complexity. Algebraic geom- 
etry told us that many geometric invariants may be computed from objects that have a 
more algebraic flavor: finite free resolutions, cohomology, Hilbert function, etc. 

There are many ways of splitting the algorithmic problems into parts, we will choose 
the following one: 

1) Provide algorithms that are easy to program, 

2) Estimate their complexity in terms of the output and/or the input, 

3) Bound the complexity of the output in terms of the input. 

As we are not at all expert in complexity theory, we will choose a measure of complexity 
that we know : Castelnuovo-Mumford regularity. It bounds the degree (in A) where most 
algebraic questions reduces to linear algebra problems. 

1. Main ingredients of two simple algorithms. 

We consider I = (fx, . . . , f s ) a homogeneous ideal in A := k[Xo, . . . , X n ] (k, a field) 
set di := deg fa and assume for simplicity that d± ^ • • • ^ d s ^ 1 and that k is infinite. 

We choose two very simple algorithms as illustrations of what we look for, there are 
more details (and other algorithms) in [Chi] for the first one and the second is based on a 
particular case of [Ch2, 5.2]. 
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The first one relies in part on the following lemma ([Chi, 20]) 

Lemma 1.1. Let gi,...,g t be a homogeneous regular sequence in A with gi = 
fi + J2j>i hijfj an d ^t J be the ideal they generate. Set a := d\ + • • • + d r — r. The 
following are equivalent, 

(1) codim(J) > t, 

(2) there exits g t+1 = ft+i+J2j>t+i h t+i,jfj such that g ± , . . . , g t+1 is an homogeneous 
regular sequence in A, 

(3) the map 

(A/ 'J). ift+ ^ fs) (A/J) a+dj 
j=t+i 

is injective. 

Steps of algorithm 1: 

- Step 1: Construct a sequence g±, . . . ,g r as in the lemma with r = codim(J), using 
(3) to determine if t = r and elementary transformations in the matrix representing this 
/c-linear map to construct gt+i as in (2) if t < r. 

- Step 2: Choose an homogeneous element h in the kernel of 

(A/ J) {fr+ ^ Js) (A/J)[d ] 

j=r+l 

of degree at most a := d\ + ■ ■ • + d r — r (there exists such an element by (3) and the 
minimal degree of such an element provides an interesting invariant of Zf. the a-invariant 
of the coordinate ring of Zj). 

- Step 3: Compute the kernel 

(A/J)^(A/J)[degh}. 

The output of the algorithm is the defining ideal of a scheme S C Zi which is purely 
of dimension dim Zj (i.e. unmixed of codimension r). If h is "general" the support of S 
is the unmixed part of Zi. The complexity is bounded by the following result (in terms of 
Grobner basis), 

Lemma 1.2. [Chi] Set b := (g u . . . ,g r , h - T de ^ h ) C A[T}. A Grobner basis 03 of 
b for the deg-rev-lex order "contains" Grobner bases for b + (h) and 1$ = b : (h). The 
maximal degree of an element in 25, for general coordinates in the Xi 's, is at most 

max{(ii + • • • + d r — r, reg(S) + deg /i}, 

except possibly for the element T 2degh . 

In characteristic zero, this bound is in fact achieved for "very general" coordinates. 

Steps of algorithm 2: Assume that I is the defining ideal of normal scheme S (i.e. 
1 = 1* and«S = Proj(^//)). 
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- Step 1: Choose two elements /, g in the Jacobian ideal of / that such that codim(/ + 
(/,<7)) = codim(J) + 2. 

- Step 2: Compute the A/ 1- module 

H\f,g;A/I) = {(x,y), fx + gy = 0}/{a(g,-f), a E A}. 

The output of the algorithm is the local cohomology module H^(A/I) ~ H (/, g; A/ 1) 
(called the Hartshorne-Rao module) . In practice choosing / and g should be easy, verifying 
the codimension condition costs quite a lot as the degrees of / and g are not that small 
when the codimension increases. 

Another strategy may be to compute first the last degree in which H^(A/I) is not 
zero (apply the same algorithm, replacing / and g by two linear forms satisfying the same 
codimension condition, the last non zero degree is the same) and then use linear algebra 
to finish the computation (in place of a Grobner basis computation, that doesn't require 
a priori bounds) . 

The main common point of these algorithmes is that they produce (at least in some 
important cases) a module that is encoding geometric informations with a complexity 
controled mainly by the complexity of the ouput. 

2. Castelnuovo-Mumford regularity 

There are many ways to define this invariant attached to a finitely generated graded 
module over a polynomial ring. Let us recall some of them in a proposition and then 
connect it to degrees of element in a Grobner basis. 

Definition. Let A be a polynomial ring over a noetherian ring k, M be a finitely 
generated A-module that is graded (for the standard grading of A), and m the ideal 
generated by the variables. For an integer i we set 

a t (M) := max{/i | W m {M)^ £ 0} 

and 

bi(M) := max{/i | Torf (M, % ^ 0} 

(with the convention max0 = — oo). 

We recall that cii(M) is indeed finite (Serre's vanishing theorem) and that the Tor 
module is equal to the Koszul homology module Hi(x;M), where x denote the set of 
variables (because K 9 (x;A) provides a free A-resolution of A/m = k as an A-module). 
Also, if k is a field and F, — > M — > is a minimal free resolution of M, then Tor^(M, fc) M := 
Hi(F, <S>a A/m) = (Fi ®a A/m)^ (maps in F, are represented by matrices with entries in 
m) is the number of minimal generators of degree fx of Fi, because A[— j] ®a A/m = k[—j] 
is concentrated in degree j. 

(See [Ei, Ch. 17] or [BH, Ch. 1] for the definition and basic facts on the Koszul 
complexes K,(z; M), K'(z; M) associated to a module M and a tuple z of elements of A. 
We will denote by Hi(z; M) and H l (z;M) their homology (resp. cohomology) modules.) 
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For simplicity, we will assume the base ring to be a field in the following, 

Proposition 2.1. Let A be a polynomial ring over a Geld k and M be a finitely 
generated A-module that is graded (for the standard grading of A). 
The following definitions are equivalent, 



(4) Let z be a finite collection of homogeneous elements of A >0 such that Mj (z)M is 
a finite dimensional vector space, 



The equivalence of (1) and (4) is easy using the standard tool for comparing two 
homological objects (spectral sequences), it implies the equivalence with (2) (taking the set 
of variables for z). Using the equivalence of (1) and (2) and the fact that a,i(M) = a^M^^) 
for any \x and % > 0, leads to the equivalence with (3). 

The definition (1), despite its apparent inaccessibility (e.g. the local co ho mo logy 
modules are not all finitely generated, unless M is of finite length) happens to be most 
tractable when one wants to estimate the regularity. As always, some familiarity with the 
object makes them very concrete and most of their apparent pathologies are not so bad 
(e.g. the graded duals of local cohomology modules are finitely generated, as they are 
isomorphic to some Ext modules). 

Definition (2) is interesting for looking at families of schemes (the study of the Hilbert 
scheme). The condition is that the maps in the minimal free resolution of over A 
have linear forms as entries. The sheaves associated to M and are the same. The 

existence of a priori bounds on the regularity is one ingredient for proving the existence 
of Hilbert schemes (see [Mu], or [EH, Ch. VI] for an introduction). 

Note that definition (4) implies in particular that if M is of dimension d and is a 
finitely generated S-module where B := k[l\, . . . , Id] and the l^s are linear forms (in other 
words, we have a Noether normalisation) then the regularity of M as a S-module is the 
same as the one as an A-module -this is also quite immediate from (1). 

One way to connect the regularity to degree of generators of a Grobner basis, is to 
study how it behaves when passing modulo a "general" linear form (see e.g. [Ei, 20.20 and 
20.21]). The key is the following lemma, where A is a polynomial ring over a field, 

Lemma 2.2. Let M be a finitely generated graded A-module and I a linear form. 
Set K := : M (I) = {m E M | Im = 0}. Then reg(M ) ^ max{reg(K), reg(M/(/)M)}, and 
if the Krull dimension of A/ Ami(K) is at most one, then 



(1) reg(M) 

(2) reg(M) 

(3) reg(M) 



maxj{ai(M) + %}, 
maxj{6i(M) - i}, 

min{^ | h{M^^) ^ b Q (M^^)+iVi} 



reg(M) := min{^ | H\z;M) > ^ l = Vz}. 



reg(M) = max{reg(K),reg(M/(/)M)}. 
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The proof is a standard diagram chasing using the local cohomology definition of the 
regularity (which has the advantage that H^(N) for i > dimiV). This lemma gives in 
particular the following, 

Proposition 2.3. Let S = Proj(A/J) be a projective surface (i.e. an unmixed 
scheme of dimension two). Assume that S fl {X n = X n -i = 0} is a zero dimensional 
scheme. Then, for the deg-rev-lex order, 

reg(I) = reg(m(J)). 

Here in(I) is the ideal generated by the leading monomials of the polynomials in / 
for a given order on the monomials. We recall that 53 is a (minimal) Grobner basis of I if 
{^ n ( f) I / £ 25} are (minimal) generators of in(I). Therefore, the maximal degree of an 
element in a minimal Grobner basis of / is bo(in(I)). The deg-rev-lex order on monomials 
is obtained by refining the degree order by the inverse of the lexicographic order. (See e.g. 
[Ei, Ch. 15] for an introduction on Grobner bases.) 

With no geometric hypotheses on Proj(A/J) one of the main early discoveries is the 
following result of Bayer and Stillman, 

Theorem 2.4. [BS] For any order and any coordinates, reg(I) ^ reg(m(/)). For the 
deg-rev-lex order, and in general coordinates, reg(J) = reg(m(/)). 

The expression "general coordinates" means that there exists a Zariski open subset 
of the linear group so that any matrix of this open subset gives rise to coordinates that 
satisfies the given property (in particular it may be that over finite fields an extension of 
the base field is needed to find good coordinates). 

Also, Diana Taylor find an explicit resolution of monomial ideals, that in particular 
proves the following, 

Proposition 2.5. If J is a monomial ideal in A, 

bi(J)^(i+l)b (J), 

so that reg(J) ^ (n + l)(b (J) - 1) + 1. 

Note that, in this form, the result is optimal: consider the case where J is generated 
by the bo-th powers of the variables. 

Mayr and Mayer, and other since then, provided examples of binomial ideals J where 
bi(J) is much bigger than bo(J) (say reg(J) 3> bo(J) c where n is the number of variables 
and c > 1 is close to V2). 

Bounds on the regularity follows from an inductive argument on the number of vari- 
ables from the following result (see [BM, the proof of 3.8]). The bounds are more or less 
of the same type as the lower bounds coming from Mayr-Meyer type examples. 

Proposition 2.6. If M is a finitely generated graded A-module, for a general linear 
form I and max{reg(M/(/)M) + 1, b (M), h(M) - 1}, 

reg(M) ^ /I + dim k H^M)^ <: /i + dim k 
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Idea of the proof. Notice that, for a general / the kernel K := ker(M — > M) is of finite 
length. Some diagram chasing gives bo(K) ^ max{6 (Af), bi(M) — 1, &2(M/(Z)M) — 1} so 
that K has no generator of degree bigger than N := max{reg(M/ (l)M) + l, bo(M), b\{M) — 
1} (in particular = implies K^ +1 = for fx ^ N). 

Now the local cohomology definition shows that reg(M / H^(M)) < reg(M/(/)M) and 

gives an exact sequence -> -> #£(M) M H H^(M)^ +1 -> for ^ > reg(M/(/)M), so 
that for fx ^ N, dim k H^(M)^ is a strictly deacreasing function of jU until it reaches 0. 

It was also remarked (and proved) by Andre Galligo that initial ideals have an interest- 
ing property in general coordinates and characteristic zero: they are stable, which means 
that if a monomial a^m is in the initial then it is also the case of xjm if j < i. Afterwards, 
Eliahou and Kervaire find a minimal free resolution for stable monomial ideals; it follows 
from the resolution that all graded Betti numbers of these ideals may be easily read from 
the minimal generators, in particular reg(J) = bo(J) for a stable monomial ideal J. 

The conditions of genericity needed for having a stable monomial ideal are not so often 
realized without performing a change of coordinates, they are more diffcult to achieve than 
the ones for having reg(J) = reg(m(/)) (for deg-rev-lex order). 

It is also important to notice that doing a generic change of coordinates have quit a 
big influence on the size of the computation, for several reasons: the coefficients get bigger, 
the polynomials became dense and the number of generators of the initial ideal increases in 
general. On the other hand it should be noted that the degrees of generators of the initial 
ideal in special coordinates may be much bigger than in general coordinates. Applying the 
following lemma to the Mayr-Mayer ideal provides such an example, 

Lemma 2.7. Let /i, . . . , f s be forms of degrees d\, . . .d s in A and set I := (/i, . . . , f s ). 

(1) If codim(J) = s, reg(J) = di H \-d s -s + l. 

(2) If codim(J) = r < s, there exists a graded complete intersection of degrees di, . . . d s 
in A[Yi, . . . , Y s - r ] such that the regularity of its initial ideal for the deg-rev-lex order 
bounds the regularity of the initial ideal of I for the deg-rev-lex order. 

The regularity of a graded ideal is always bounded in terms of the Hilbert function of 
the ideal, as all the graded Betti numbers are bounded above by the ones of the lex-segment 
ideal that only depends on the Hilbert function (see [Bi] and [Hu] for characteristic zero 
case, [Pa] for the general case, and [CGP] for a short argument in characteristic zero). If / is 
saturated, the regularity of the lex-segment ideal only depends on the Hilbert polynomial. 
The regularity of the lex-segment ideal may be computed ([CM, 1.3 and 2.3]) and leads 
for example to the following bound that is at least as bad as expected... 

Corollary 2.8. [CM] Let I C A be an ideal generated by polynomials of degrees 
di, . . . ,d s and let r be the codimension of I. For any admissible monomial order and any 
coordinates, 

reg(m(/)) < 1 + [di • • • d s f l ~ r 
ifr^n and reg(m(/)) ^ d\ + • — h d n+ \ — n if r = n + 1. 
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3. Bounds on Castelnuovo-Mumford regularity 

There is a famous conjecture that suggests the following bound for reduced and irre- 
ducible schemes: 

ConjecturefEisenbud and Goto]. If S C P n is a non degenerate reduced and irre- 
ducible scheme, 

reg(«S) ^ degS — codiuuS. 

(Non degenerate means S <£. H for any hyperplane H.) 

We recall that if S := Proj(A/J), reg(<S) := reg(A/P) = reg(P) - 1. 

This result was known for curves when the conjecture was made. It was first es- 
tablished for smooth curves by Castelnuovo [Ca], and the for reduced curves with no 
degenerate component by Gruson, Lazarsfeld and Peskine (over a perfect field) in [GLP]. 
There is some evidence that this may be true at least for smooth schemes in characteristic 
zero: it is true for smooth surfaces (Pinkham and Lazarsfeld) and (up to adding small 
constants) in dimension at most six, by the work of several people including Lazarsfeld, 
Ran and Kwack. 

In any dimension, it was prove by Mumford ([BM]) that in characteristic zero a smooth 
scheme S satisfies, 

reg(S) < (dim S + l)(deg«S - 1). 

In positive characteristic, it follows from theorems that we will mention below that 
one has reg(<S) ^ (dim5 + l) 2 (deg£ — 1), and there are also quite reasonable results for 
schemes with isolated singularities. 

We now turn to bounds depending on the degrees of generators. As we mentioned in 
the preceeding paragraph, there is no reasonable bound on the regularity without imposing 
geometric conditions. 

Let I = (fi, . . . , f s ) be a homogeneous ideal in A, where fi is a form of degree d^. 
We will assume that d\ ^ d^ ^ • • • ^ d s ^ 1. Let Zj C P n be the scheme defined 
by I and r be the codimension of I in A, which is also the one of Zj as a subscheme 
of P n . Let S be the top dimensional part of Zi and Y the residual of S in Zj. In 
algebraic terms, is is the intersection of the primary components of I of codimension r, 
and I Y := (I : I s ) = {/ e A | fl s C I Y }. 

The first striking result on regularity in these terms is due to Bertram, Ein and 
Lazarsfeld: 

Theorem 3.1. [BEL] If Zj = S is smooth of characteristic zero, 

reg(S) ^ d\ H h d r — r, 

with equality if and only if S is a complete intersection of degrees d±, . . . , d r . 
This linear bound was generalized in [CU], 

Theorem 3.2. Assume that S have at most a one dimensional singular locus an is 
locally a complete intersection outside finitely many points. If the residual Y have at most 
isolated singularities and k is of characteristic zero, 

reg(S) < d 1 H h d r - r, 
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with equality if and only if Zj = S is a complete intersection of degrees di, . . . , d r . 

Note that the defining ideal of S may be computed in low degree by Algorithm 1, even 
if the regularity of / is much bigger. They rely on liaison theory and use either Kodaira's 
vanishing theorem or a result of Karen Smith that enables an induction on the dimension. 

More recently, we showed several other bounds. They essentially improve the ones of 
[CU] in positive characteristic, and provide the following result: 

Theorem 3.3. [Ch2, 4.4] Assume that Z is an isolated component of Zj that doesn't 
meet the other components, and that Zj is smooth at all but a finite number of points of 
Z. Then, 

reg(Z) ^ (dimZ + l)(d\ H h d r - r - 1) + 1. 

This generalizes the result of [CP] that treats the case where dim Z = 0. The proof 
relies on [CP] and a result of Hochster and Huneke, which implies that the phantom 
homology (which is, roughly speaking, the one that vanishes in the Cohen-Macaulay case) 
is uniformly killed by the Jacobian ideal. The result then follows by cutting the scheme 
Z by a sequence of parameters in the Jacobian ideal and using some homological algebra 
to exploite this uniform vanishing. The connection between annihilators and vanishig 
was already remarked and used to study the so-called £-Buchsbaum schemes by Miyazaki, 
Nagel, Schenzel and Vogel ([Mi], [NS1], [NS2] and [MV]). 

Let us also point out the following remark that formalizes the fact that bounding the 
regularity in a geometric context is as difficult as bounding the degree where the Hilbert 
function becomes a polynomial, or bounding the degree where every global section is the 
restriction of a polynomial. 

Remark 3.4. [CM 2.5] Let V be a property of embedded projective schemes and N(X) 
a numerical invariant attached to such a scheme X. Assume that if X C P n satisfies V and 
H is a general hyperplane, then X n H C H ~ P n _i satisfies V and N(X n H) ^ N{X). 
We denote by Ix ^ R the defining ideal of X and by Hx the Hilbert function of R/Ix- 
Then the following are equivalent, 

(i) If X satisfies V, veg(R/I x ) < N(X). 

(ii) If X satisfies V, reg(H x ) < N(X) - 1. 

(iii) If X satisfies V, {R/Ix)» = H°(X, O x {p)) for /i ^ N(X), 

where reg(H x ) is the last degree where H x differs from the Hilbert polynomial Px- 

Examples for property V are: X satisfies Sk, X is smooth in codimension £, X is 
irreducible, X is equidimensional, or any conjonction of some of these properties. For N(X) 
one may choose the degree of X , or the degree of X minus the embedding codimension of 
X if X is irreducible and reduced, or the minimum over the sets of equations defining X 
of the maximal degree of these equations. 

Another important point is to notice that even if we are not able to bound the reg- 
ularity in many cases, a big part of the information is sometimes available in an indirect 
way. For example, if the top dimensional component S of Zj have at most isolated sin- 
gularities the canonical module us of S have a small regularity (at least in characteristic 
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zero, thanks to Kodaira's vanishing theorem) and is easily computable (as the kernel of 
the map in Step 2 of Algorithm 1). From u>s we may compute the Hilbert polynomial of 
S or its cohomology modules using Serre duality (at least if S is Cohen-Macaulay), or test 
if an element is in 1$. This in turn gives a way to check if reg(5) ^ N by linear algebra 
computation in degree at most N plus a linear function of the degrees of generators (in 
general coordinates, the criterion of [BS] gives such a test). 
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